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Abstract 

The phase diagram of a model describing doped CuGeC>3 is derived. The model emphasizes the 
role of local moments released by the impurities and randomly distributed inside the gaped singlet 
background. The phase diagram is investigated by two methods: (i) in a mean field treatment of 
the interchain coupling and (ii) in a real space decimation procedure in a two dimensional model of 
randomly distributed moments. Both methods lead to similar results, in a qualitative agreement 
with experiments. In particular, a transition to an inhomogeneous Neel phase is obtained for 
arbitrary small doping. From the decimation procedure, we interpret this phase at very low 
doping as a Griffith antiferromagnet. Namely, it does not have a true long range order down to 
zero temperature. Nonetheless, large magnetically ordered clusters appear already at relatively 
high temperatures. This demonstrates the role of disorder in the theoretical description of doping 
in CuGe03. A detailed comparison with other approaches is also given. 
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1 Introduction 



Spin-Peierls systems have been the subject of a renewed interest because of the discovery of several 
inorganic quasi-one dimensional (Q1D) spin-Peierls compounds. A very well studied example is 
CuGe03. This Q1D compound made of weakly coupled CuC>2 chains has a spin-Peierls transition at 
Tgp ~ 1AK O, as revealed by diffraction measurements, showing a dimerization below Tgp, and by 
susceptibility data, which indicate the opening of a spin gap at Tsp. 

CuGe03 can be doped in a controlled fashion, giving the unique opportunity to study the role of 
impurities in a spin-Peierls system. Doping is done in two different ways. One possibility is to directly 
affect the spin degrees of freedom within each Cu02 chain by substituting some of the Cu atoms 
(whose spin is S = 1/2) with magnetic (Ni ||], which has S = 1, or Co p), S = 3/2) or non-magnetic 
(Zn |j| ||], Mg |7|) impurities. Another way is by substituting Ge with Si ||, which is expected to 
modify the strength of the exchange couplings among the magnetic ions. For all these compounds, 
even a small amount of doping strongly enhances the magnetic fluctuations, leading in most cases 
to a lower temperature Neel phase coexisting with the Peierls's distorted phase. For instance, in 
0.3% Si-doped compounds, free spins are released which contribute with a Curie-like behavior below 
Tsp — 12. 5K (which is slightly reduced with respect to the undoped sample), and finally get frozen 
below the Neel temperature T\r — 0.95/^ ||. The same behavior has been reported in a 2% Zn-doped 
compound, where magnetization measurements can be interpreted as if approximately one spin-1/2 
per Zn atom is released IfTcfl . In all cases but Co-doped samples [11], the low doping x < 1% phase 
diagram seems to have universal features [12]: T$p decreases almost linearly with x, while no critical 
concentration for the appearance of the Neel phase has till now been revealed down to the lowest 
accessible impurity concentrations (e.g. 0.12% Zn for which T/v — 0.025^ [|l4|]). 

At larger doping, the situation is still not fully established. However, most of the data indicate 
that over some critical x c , the spin-Peierls long-range order disappears, although short-range corre- 
lations seem to persist. For Mg-doped samples, it is claimed that at x c ~ 2.3% a first order phase 
transition takes place which is signaled by a jump of |7j. 

Manabe et al. [14] have reported an antiferromagnetic (AF) ordering with Zn concentrations as 
low as 0.11%. Their analysis of the variation of the Neel temperature versus doping concentration has 
shown the absence of a critical doping concentration for AF ordering at zero temperature {x c = 0). 
The goal of the present article is to understand this onset of long range AF order at such low doping. 
This experimental observation questions the validity of existing theories in a much more drastic way 
than the coexistence of AF order and dimerization. 

Several distinct theories describe this coexistence of dimerization and antiferromagnetism. Fukuyama 
et al. |3q] first proposed a model showing that the two orders can coexist, and another model was 
proposed by Mostovoy et al. f27j]. Two of the present authors have proposed another model for the 



coexistence of AF ordering and dimerization [16, 17]. All these theories successfully describe the 
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coexistence of dimerization and AF order upon doping, and they also show that the effect of doping 
amounts to generate states inside the Spin-Peierls gap, in agreement with susceptibility experiments 
(for example A similar behavior is reported in spin-ladder systems, which, in spite of having 

a large spin gap, undergo a Neel transition upon very small doping. For instance SrCu203, with a 
spin gap of ~ 400K, has a finite T/v — 2>K already at 1% Zn doping |l~5|| . In these ladder compounds 
the low energy spin excitations induced by doping have also been revealed by a finite specific heat 
coefficient [jpjfl . 

In their article, Manabe et al. pointed out that one energy scale emerging from our model 
(denoted by T*, see sections |2.2| and || for its definition) was many orders of magnitude lower than 
the Neel temperature (for instance at x = 0.1%, T* ~ 4 x 10 _42 K whereas T/v — 20mK has been 
measured p^l ). On the basis of the analysis of AF fluctuations in the ID model, we already suggested 
[16, 17] that this apparent inconsistence can be avoided if one makes a proper treatment of disorder 



averaging. In the present article, we push further this analysis by estimating the Neel temperature 
in the presence of interchain coupling, which turns the quasi-long-range AF fluctuations of the ID 
model into a longer range order. We prove by two different methods (a mean field treatment of the 
interchain coupling, and a decimation procedure of the 2D disordered problem) that the resulting 
Neel temperature has the right order of magnitude. Moreover, under both treatments, we show that 
there is no critical concentration for the establishment of AF order. As it was shown by Mostovoy et 
al., this is not the case in their model, and we believe that this is not the case either in the approach 



by Fukuyama et al., further extended by Yoshioka and Suzumura 3q| . 

The paper is organized as follows. We introduce the model in section ||, and discuss the experi- 
ments by Manabe et al. |l4j] . The problem of interchain coupling is first solved in section [|| at the 
mean field level. Section ^| is devoted to the decimation approach that goes beyond the mean field 
treatment of interchain interactions, and is more relevant for the extremely low doping situation. The 
picture of Griffith- like antiferromagnetism emerges from this treatment: large AF correlated clusters 
are formed when the temperature is decreased below the spin-Peierls temperature, without infinite 
range ordering. A summary of our results is presented in section |5| together with a comparison with 
the approaches by Mostovoy et al. p7|], Fukuyama et al. [35], and Yoshioka and Suzumura [36]. We 



argue that these approaches would fail to properly describe the low doping behavior, whereas our 
description leads to AF ordering at infinitesimal doping. Since the analysis of the work by Fukuyama 
et al. is more technically involved, we have kept the qualitative features for the concluding section 
while the technical aspects are left for the appendices. 

2 The Model 
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2.1 Generation of spin- 1/2 moments 

We start by briefly describing our model. Let us consider the Hamiltonian of a dimerized Heisenberg 
model. Given the large anisotropy in the couplings (for instance f[§| l) 

\J a \ ~ J c /100 < J b ~ J c /10 < J c ~ 10.6meV, (1) 

where c is the chain axis, and a and 6 the two perpendicular directions, the Neel temperature is 
mainly determined by J c and Jb. We therefore consider the Hamiltonian of a dimerized Heisenberg 
model in two dimensions, 

1 + J± S n ,i • S n +l,i, (2) 
n,i n,i 

where J = J c , J± = Jb, i labels the sites along the chain and n is the chain index. 5 is the 
dimerization parameter which is staggered along the b axis. In reality, frustration has to be added to 



get a reasonable agreement with the spin susceptibility data above T$p [21|. This can be accomplished 
by introducing for instance a (quite large) next nearest neighbor exchange J' ~ 0.23 — 0.36J, which 
can result both from next nearest neighbor Cu-Cu superexchange and from the lattice quantum 
fluctuations. However, this additional complication is neglected in our present analysis. 

When 5 > 0, each spin at an even site is coupled by a strong bond J(l + 5) to the spin at its 
left, and by a weak bond J(l — 5) to the spin at its right, both on the same chain. A finite gap 
between the singlet ground state and the lowest triplet excitation arises as soon as 5 ^ and J± 
is less than a critical value which depends on 5. This should be the case of the pure CuGeOs. We 
start by assuming that the most important effect of doping is that each impurity releases a free 



spin out of the dimerized singlet background. This behavior was anticipated in |22| on the basis 



of exact diagonalization of small-size systems, and in Ref. [16] by more qualitative arguments (see 
Fig. ||). Further numerical works have confirmed this picture |23|, 24], and have shown its robustness 
when more realistic assumptions are made, such as quantum phonons or coupling to a 3D lattice 
dimerization [25], or next-nearest-neighbor frustration [26|. The generation of spin-1/2 objects is 



most clearly seen for Zn or Mg doping and in the limit J± = 0. As a first approximation, the non 
magnetic ion cuts the chain into two segments, e.g. one starting with a weak bond, the other ending 
with a strong bond. In the former segment, a free spin-1/2 moment will appear localized at the 
boundary. Even if we add a weak link across the non-magnetic ion, which couples the two end-spins 
of each segment, this free spin-1/2 moment still exists. An even more transparent explanation is 
obtained by considering a ring of even number of sites. Since 5^0, the ground state is a singlet with 
a gap to the lowest S = 1 excitations. If a spin is removed, and a weak link across the empty site 
is added, the effective spin chain will now contain an odd number of spins 1/2. Hence the ground 
state will have S = 1/2, in spite of a finite gap to higher spin states. Moreover, since translational 
symmetry is broken, this S = 1/2 moment will be mostly localized around the weak link. Similarly, 
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Figure 1: Mapping of a segment with a non magnetic site to a squeezed segment, without this site. 
The double line represents a strong bond in the dimerized chain, i.e. J(l + 5), the single solid line a 
weak bond J(l — S), while the dotted line is the weak link across the non magnetic site. Going from 
(b) to (c) we have assumed the weak link equal to the weak bond. Therefore the squeezed segment 
in (c) contains a domain wall, i.e. two consecutive weak links. 



even if one spin 1/2 is substituted by a spin 1, modeling Ni doping, a free spin 1/2 moment is formed 
around the impurity. We also believe that the same situation holds with a finite Jj_. 

The case of Si-doping is less clear. Ge is supposed to influence the superexchange between the 
Cu's fll|, increasing the antiferromagnetic superexchange with respect to the ferromagnetic indirect 
exchange, finally leading to a net J > 0. The smaller bond length of Si-0 has likely the effect 
of reducing the antiferromagnetic superexchange by diminishing the angle Cu-O-Cu, thus inducing 
defects in the dimerized structure. In principle, a large reduction of that angle might even change 
the sign of the Cu-Cu exchange constant, leading to the formation of free spin-1 moments. Although 
we believe this is likely to occur, here we will limit our analysis to the simpler case of non magnetic 
ions substituting Cu. 

We therefore take for granted the appearance of spin-1/2 moments upon doping, in the line 



of several numerical works (see e.g. Refs. [22, 23, 24, 25, |26|) and pursue the next step which 



is to build up the phase diagram of doped CuGeOs from these excitations. This cannot be done 
without including the effect of interactions between those spin-1/2 moments, which are central to 
our description. 

As we argued, in such cases free spin-1/2 moments form around each impurity, showing up as low 
energy states inside the spin-Peierls gap, which, at very low doping, is almost unchanged. Therefore, 
we model the doped compounds as a two components system, where spins 1/2 are diluted inside 
a singlet background, characterized by a fixed gap between the singlet ground state and the first 
excited triplets. The virtual polarization of the singlet background provides an antiferromagnetic 
coupling between these spins. Since the dimerized state is gaped, this coupling decays exponentially 
with the distance between two spins over a length of the order of the spin-Peierls correlation length. 
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For instance, let us consider two spins at a distance r = (n x c, n y b), n x and n y being positive integers, 
and c and b the lattice constants in the chain and perpendicular to the chain directions. Then, the 
exchange between these two impurity-released spins can be taken of the form 

/ 

J(f)~ -(-!)"«+"» A exp - 

V 



where A is the spin gap, £, x = £,sp — 9 — 13 the correlation length in lattice constant units along 
the c-axis, while £ y ~ £spJ±/J — 0.1£sp, the correlation length along the b-axis. The oscillating 
factor in front takes into account the fact that the singlet background polarizes mostly in a staggered 
way, leading to a ferromagnetic exchange if the two spins are on the same sublattice and to an 
antiferromagnetic one in the opposite case. Since the impurities are located at random, replacing a 
fraction x of the magnetic sites, Eq.(^) defines a two dimensional random Heisenberg model, which 
should properly describe the excitations at energies below the spin-Peierls gap. 

The above two-component model looses its validity at larger doping, where the impurities start to 
affect the singlet background in a non negligible way. Nevertheless, we believe that at low doping it 
does capture the essential physics of the real system, giving a good description of the the Neel tran- 
sition. Below some temperature, these weakly coupled spins will tend to order antiferromagnetically. 
Moreover, they will also induce a staggered polarization on the singlet background, which does have 
a finite staggered susceptibility, in spite of an exponentially decaying uniform susceptibility. Accord- 
ing to this scenario, the Neel phase would consist of spins antiferromagnetically ordered diluted in a 
partially polarized singlet background. The main problem is to get an estimate of the Neel transition 
temperature, which allows a direct comparison of the present model with the experimental phase 
diagram of Cui_ x Zn x Ge03 fi"4j| , where it has been shown that zero temperature AF ordering occurs 
at an infinitesimal Zn concentration. 

2.2 The typical energy scale T* and the susceptibility experiments by Manabe et 



al. [P 



We first proceed by considering the extremely low doping susceptibility experiments by Manabe et 
al. [14], which, to our knowledge are the only available with Zn concentrations as low as 1.12 x 1CP 3 . 



As pointed out in Ref. [14], the Neel temperatures (see table [l]) are indeed much larger than the 
typical energy scale T*, given by the Neel temperature of a regular array of equally spaced spin-1/2 
moments with an exchange given by Eq. (Q), which we may estimate by considering, instead of a 
random distribution, a regular distribution of the spins over the 2D lattice with a concentration x. 
The coupling along the c and 6-axes are then J x = A exp(— l/£ x ) and J y = A exp(— l/Cy), respectively, 
with I = \j \fx the average distance between the spins. Hence, at low doping, the effective Heisenberg 
model describing the spins has a huge space anisotropy J x /J y ~ exp(— l/£ y ) <C 1, implying that the 
scale controlling the Neel temperature is T* ~ J y ~ 2 x 10 _14 K if x = 0.1%, which is far too small 
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Figure 2: Fit of the susceptibility experiments by Manabe et al. [14| to the Curie form (||) (solid 
lines) and to the form @ (dashed lines). For clarity, we have mutiplied the susceptibility by the 
concentration, so that all the concentration data are at the same scale at low temperature. Since 
the data provided to us by Manabe et al. were very dense in some temperature regions, we have 
"compressed" the data by averaging them over given bins. The parameters of the fits are shown on 
Table |. 

compared with the experimental values of the AF transition temperatures. The typical energy scale T* 
can also be evaluated in the single chain version of our model, as recalled in section ||. Experimentally, 
Manabe et al. |l4j measure 0.027K< TV < 0.7K in CuGeOs, for concentration impurities ranging 
from 0.11% and 0.49% (see Table [l]). A more accurate estimate of the characteristic exchange 
interaction can be made by fitting their data in the range T/v < T < T$p. 

We first decompose the susceptibility in the whole range between and Tgp into two contribu- 
tions: 



X(T) 



a 



imp 



exp 



Ji 



imp 



+ 



Csp 



T exp.-- 



A 



(4) 



T \ T 

the first term corresponding to the low temperature contribution of Zn impurities, and the second 
term describing the higher temperature behavior of the susceptibility close to the spin-Peierls tem- 
perature. The energy scale involved in this contribution is the spin-Peierls gap A ~ 44.7 K. The fits 
for various impurity concentrations are shown on Fig. ||, and are in a quite satisfactory agreement 
with experiments. 

The three parameters Ci mp , Ji mp , Csp are shown on Table [|. The energy scale Jj mp involved 
in the fit of the impurity contribution sets the characteristic energy scale involved in the physics 
of impurity-released magnetic moments. As it can be seen on Table |], Ji mp is finite, which shows 
that the Zn-induced magnetic moments are interacting, consistent with our picture. Moreover, this 
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Table 1: Parameters of the best fits of the susceptibility data by Manabe et al. to the forms (0) 
and (||) shown on Fig. |2|. The Neel temperature is also given, and is of the same order of magnitude 
as J; mp and 9. The values of << J >> (Eq. (||)) and T* (Eq. (||)) are also given. 



energy scale is not of the order of the typical exchange T* . 



Alternatively, the susceptibility data in [14] can be fitted with the Curie- Weiss expression 



6> + T T r v T 

a fit proposed in Ref. ||. As it is visible on Fig. ^, the form of the susceptibility also leads to a 
good description of the susceptibility for temperatures ranging between the Neel temperature and 
the spin-Peierls temperature. The parameters of the fits are shown on Table |l], and the energy scale 
9 involved in the fit (||) is again orders of magnitude larger than the typical exchange T*. Both 
expressions (0) and (||) have the same high temperature behavior in the limit T 3> J, 9, and, in this 
regime, 9 = J. Once the two fits have been performed independently, we can check that these energy 
scales are indeed the same, as it is visible on table ED, and comparable to the ordering temperature 
T N . 

The reason why the Neel temperature is not of the order of T#, but many orders of magnitude 
higher, is that, as a result of disorder in the spatial distribution of magnetic moments, clusters exist 
in the system with spins which are close to each other, being therefore strongly AF correlated. Once 
the interchain couplings will be included (sections ^ and |j), this will result in a Neel temperature 
many orders of magnitude larger than T* . 

3 Mean field theory of AF ordering 

In this section we first recall (section |3.1| ) the results obtained in ID [ Iq , 17], discuss the energy scale 



T* in the ID model (section |3.2|), and present in section 3.3 the extension to include the effect of the 



interchain coupling Jj_ . We take in this section as well as in the forthcoming section |] a spin-Peierls 



gap A = 44. 7K [18], assumed to be independent on temperature. The spin-Peierls correlation length 
is Csp = 9 M. 
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3.1 Physics of the ID model 



As shown in [16, [17]], a single chain does develop quasi- long range magnetic correlations at zero 
temperature, which we believe is the key point for understanding the establishment of AF ordering at 
arbitrarily low doping. More precisely, the correlation length £(T) is found as the root of f(y, IY)+1 = 
0, with y = -l/£, T T = hi (A/2 1 ), and 

fo(y) - atan(ar) 



f(y,T) = a 



/o(y)tan(«r) -fa' 



where a = £,sp\/—y(y + 2xJ, and fo(y) = £sp(u + %)■ The same solution holds in the XX limit, 
except that the correlation length is then found as the root of the equation f(y,TT) + 2 = 0. The 
correlation length is compared on Fig. || to that of an ordered array of magnetic moments separated 
by a distance of 1/x and interacting with an exchange A exp (— 1/x^sp)- 
Finally, we will use latter the local staggered susceptibility Xs(T) 

Xs{T)=X { X ln(A)T)ll/t SP ) W (6) 

3.2 Energy scale T* 

The staggered spin-spin correlation function 



(S>(r)S*(p)) ~ (-1) 



shows quasi-long-range AF fluctuations at low temperature, with the limiting behavior of the corre- 
lation length 

^ X ^SP i„2 (T 



OT^^ln^-j (7) 

valid at extremely low temperature T < 2 1 *, in which case the system has crossed-over to the random 
singlet fixed point |l7], . The distribution of the exchange constants is 

This distribution P(J) is such that the typical exchange constant @, when x^sp 1, is much 
smaller that the average exchange 

« J(r) »=A - X$5 f . (8) 

The energy scale is defined by 

~ e «iog[J r (r)]» _ / ^ e -i/fep) ) (9) 
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Figure 3: Variations of the correlation length in the XX limit for x = 0.3% doping versus temperature, 
with £sp = 9 an d A = 44.7 K. (a) shows the variations of the correlation length of the disordered 
system and (b) the correlation length of the typical realization of disorder, namely, an array of 
equally spaced magnetic moments. The limiting value of the correlation length when T — > T$p is 
l/(n + 2/£gp), the factor of 2 originating from the XX approximation. 



and is equal to the typical exchange constant among the spins. The symbol << ... >> denotes a 
disorder average [ 32 1 . The definition of T* in Eq. (0) in the ID model is equivalent to the one we 



have introduced in section 2.2 on the basis of a 2D modeling. Above this temperature, at which the 



correlation length is of the order of 1/x (the average distance between the impurities), £(T) decreases 
with increasing temperature, approaching £$P when T —> A. As shown on Fig. ||, at intermediate 
temperatures T* < T < Tsp the AF correlation length is increased by disorder. Above T*, the 
disordered chain takes advantage of rare fluctuations in which many spins get closer than the average 
distance thus forming AF correlated clusters due to the exponential dependence of the exchange 
(||) upon the distance. 

The experimental values of the Neel temperature (see Table [l]) are clearly of the order of the 
average exchange << J >>, and not of the order of the typical exchange T*. The same is true for 
the characteristic exchange Ji mp in (|j) and the (negative) Curie temperature 9 in (||). The orders of 
magnitude of Ji mp and 9 obtained in section |2.2| are in fact a proof of the relevance of disorder in the 
physics of the system. 
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Figure 4: Comparison between the experimental phase diagram by Manabe et al. [jl4[ (O), our 
mean-field treatment (+) and the decimation results (L = 40: □; L = 81: x). The fit to the 
experimental data is ^4exp (—B/x), with A = 2.3 K and B = 5.7 x 10 -3 . The mean field results (+) 
have been obtained with £sp = 9 c and Jj_ = 12 K. The decimation results have been calculated 
with sizes L = 40 (□) and L = 81 (x), and correspond to an existence probability of 1/2. As it can 
be seen on this figure, both estimations of the Neel temperature at low doping are much larger than 
T* = Aexp (— 1/x^sp) (— 10~ 47 K if x = 0.1%). The experimental results as well as the results from 
the mean field treatment and the decimation procedure show no critical concentration. The mean 
field and decimation results have been obtained by using the known parameters (|l]) for CuGe03 as 
an input of the calculation, without adjusting any parameter. 



3.3 Mean field treatment 

We now include the effect of an interchain coupling J± ~ 12K [|l8| under the form of a mean field 
treatment based on the ID model proposed and analyzed in Refs. [16, 17]. As recalled in section 
|3.1| , the ID version of our model does develop quasi-long-range AF fluctuations, which is the reason 
why AF order will be easily generated once an interchain coupling J±_ is taken into account. When 
the interchain exchange J± is switched on, a true Neel long range order might develop. The question 
whether the Neel order is a true long range order or a short range order will be discussed in more 
details in the light of the decimation calculation in section ||. We look in the present section at 
the mean field response to a uniform staggered field, and therefore the Neel order is infinite ranged, 
which is not the case when the interchain coupling is included beyhond the mean field treatement. 
To get an estimate of the Neel temperature, we use the following mean- field like equation for Tjy: 



J±Xs(T N )^T J 



(10) 
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where \s (T) is the local staggered magnetic susceptibility (||) , and £ (T) the correlation length drawn 
on Fig. ||[ The Neel temperature calculated via Eq.([lO|) is plotted on Fig. |I| as a function of the 
impurity concentration x for different values of Jj_ . It is quite clear from the present calculation that 
the Neel temperature T/v is much larger that T*. This originates from the fact that large magnetically 
ordered clusters form already much above T*. This fact, together with the finite polarizability of the 
singlet dimerized background, explains the surprisingly high Neel temperatures T/v. 

We can even obtain a rougher estimate of the Neel temperature. Since T* <C T/v < T$p, we can, 
in a first approximation, take the correlation length £(T) equal to the spin-Peierls correlation length 
£sp, and further assume that the susceptibility amounts to a paramagnetic contribution of the free 
magnetic moments: Xs{T) — x/T. Using Eq.(|i~0|) we obtain the Neel temperature 

T N ~ J±x£sp, (11) 

linear in the impurity concentration. Taking J± ~ 12 K, and £sp — 9c and x = 0.01 we get 
T/v = 1.08 K. At the same concentration, Martin et al. measure T/v — 2K ||, of the same order of 
magnitude as our estimate. 

Both Eqs. ( |T0| ) and ( JTT| ) are approximate expressions, which might work at intermediate impurity 
concentrations, but whose validity is doubtful at very low doping x. Anyhow, those results would 
imply that no critical concentration is needed to get a Neel phase, since a single chain does develop 
quasi-long range correlations for any non zero doping. Moreover, the Neel temperature obtained in 
this mean field treatment is obviously not of the order of T* , but of the same order of magnitude than 
the experimental value. At this point, we judge the mean field treatment to compare satisfactory to 
the experimental data, given the fact that our model is certainly schematic. These results will be 
further confirmed by the decimation analysis in section |j. 

4 Low doping physics: decimation in 2D 

We now present a decimation scheme in 2D which allows to describe the physics at extremely low 
doping. We have already shown in section |3| that our model leads to correct orders of magnitude of 
the Neel temperature. We do not believe this mean field theory to be valid at very low doping, since 
lowering the doping should enhance inhomogeneities, and we then do not expect the physics to be 
captured by a homogeneous staggered field involved in the mean-field Stoner criterion (|i~0|). This is 
why we now consider the full disordered problem in a fixed dimerized background. We take a square 
lattice of size L x L, representing the initial Cu sites of CuGeC>3, which, without doping, are paired 
into a spin-Peierls dimerized phase. Some of these lattice sites will be occupied by non magnetic 
impurities with a probability x, leaving in their vicinity unpaired spin 1/2 moments. The exchange 
coupling between two such spins at positions (xi,yi) and {x2,y2) is ffi-2 = J\-2^i • S2, where we 
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A being the spin-Peierls gap A = 44.7 K, and £ x ~ 9c and £ y ~ 0.1 £ x the correlation lengths in 
the c and 6 directions, respectively. This model has no adjustable parameter. The exchange Ji_2 
is negative if the two sites belong to the same sublattice, and positive otherwise, which takes into 
account the finite response of the singlet background to a staggered field. This defines a disordered 
spin model in 2D with long range interactions, in the sense that any spin has an exchange with 
any other, which however decays exponentially with their separation. The decimation procedure we 
use here is similar to the one used by Bhatt and Lee |[3| in the study of phosphorus doped silicon. 
However, we decimate strong bonds without renormalizing the exchanges, which is just sufficient to 
get the order of magnitude of the ordering temperature. 

By construction, our model has bonds which are never frustrated and are perfectly compatible 
with a long range AF order, being ferro on the same sublattice and antiferro on different sublattices. 
Therefore we may not get the same low temperature properties found by Bhatt and Lee in their 
work on Si:P. The reason of the difference is that our spins are diluted into a matrix which is easily 
AF polarizable. A similar physics does not exist in Si:P. 

4.1 Decimation in 2D 

In one dimension, a RG transformation of a random AF spin chain consists in picking-up the strongest 
exchange of energy E, decimate the two corresponding spins, and renormalize the exchange by 
projecting onto the singlet states of the decimated bonds. This leads to broad renormalized exchange 
distributions and to asymptotically exact low-energy physics |20| ]. 

Similarly to the study on Si:P by Bhatt and Lee |3^1 , we carry out the same procedure in 2D, 
with the difference that we do not define a renormalized exchange but just eliminate the strong 
bonds. This approximation is expected to lead to qualitatively correct results as far as the estimate 
of the ordering temperature is concerned. Considering the 2D disordered problem defined by the 
Hamiltonian (|l2|) , we decimate the bonds of strength E < J < A. Assuming the initial bonds 
were colored in white, we color in black these strongly coupled bonds. As the energy E decreases, 
more and more white bonds will be turned into black, until the black bonds percolate throughout 
the system at an energy Ep. We use Ep as an estimate of the Neel temperature Tjy, even though, 
strictly speaking, Ep is an upper bound of the Neel temperature. In practice, we use the following 
techniques to solve this problem numerically, (i) The decimation is carried out by quick sorting the 
exchanges, (ii) A linked list is built at each energy that encodes the graph of existing exchanges, 
(iii) A Hoshen-Kopelman relabeling algorithm |34|] is used for cluster labeling. We have tested our 
program with standard 2D site and bond percolation, with excellent results. We now present our 
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results, by first showing in section 4^ the failure to find a percolating cluster in the infinite volume 
limit. We next consider in section O ordering in a finite coherence volume. 



4.2 Absence of self similar ordering 

In this section, we fix the impurity concentration x and look for percolation in an initial lattice of 
size Lx L, with increasing sizes L. Given the anisotropy of the problem (£ x = 9 = 10^), percolation 
occurs along the x direction of the strongest exchanges. It may happen that the left-most and right- 
most columns contain no impurity at low concentration, because of the finite probability Lx to have 
an impurity present in a given column. We are thus lead to consider as percolating a cluster that 
percolates between the left-most and right-most columns that contain at least one impurity. 

In order to determine the percolation probability, we consider two (equivalent) quantities. First 
the existence probability El, i.e. the probability to find a percolating cluster in a system with a finite 
size L. We also consider the Binder cumulant Bl = (S 2 )l/{S) l , with S the number of impurities in a 
percolating cluster. We analyze in this section the finite size scaling behavior of these two quantities. 
More precisely, if a true percolation transition would occur (as it is the case in the 2D bond and 
site percolation problem), the crossing points either of El(T) or of Bl{T) for increasing sizes would 
accumulate at the transition temperature. We have plotted on Figs. || and || the existence probability 
and Binder cumulant for increasing sizes at a fixed impurity concentration x = 3%. As it is visible 
on these figures, there is no true percolation transition in the sense that while the crossing point 
temperatures are higher and higher as the sizes increase, the existence probability decreases to zero. 
We interpret this result as the lack of true long range order in an infinite system. Therefore, we 
believe that AF ordering is of a Griffith type. 

In fact, this finding is not unexpected. At a fixed energy T, we can interpret the diluted spins as 
particles with a core of area a approximately given by 



The core is indeed elliptic, with a ratio between the axes a/b = £,x/£,y Hence, the ratio between the 
area occupied by these particles and the surface sample is 



Ignoring the complication due to the anisotropy, and assuming that the cores are impenetrable, 
a percolating cluster can only appear if T < Tp where n e ff(Tp) = n c , being n c = 0.59 the site 
percolation threshold of a 2D square lattice. However, the core of each particle is not impenetrable. 
For this reason the existence probability does not approach one for increasing sizes even below Tp. 
Indeed, the existence probability E(T) for L — » oo does not tend to a step function as for standard 
percolation, but to a smooth curve as it is visible from the finite size study shown on Fig. || 




(13) 



n e ff(T) = xcj(T). 



(14) 
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Figure 5: Variations of the existence probability E(T) versus temperature T at a fixed impurity 
concentration x = 0.03 and for increasing sizes L = 40 (O), L = 57 (+), L = 81 (□), L = 115 (x). 
As it is visible on this figure, there is no self-similar fixed point. 
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Figure 6: Variations of the Binder cumulant B(T) versus temperature T at a fixed impurity concen- 
tration x = 0.03 and for increasing sizes L = 40 (O), L = 57 (+), L = 81 (□), L = 115 (x). As it is 
visible on this figure, there is no self-similar fixed point, in agreement with the existence probability 
result shown on Fig. ||. 
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4.3 Low doping phase diagram 

The previous analysis suggests that a true Neel long range order is not established at very low 
doping, where the singlet background is assumed only to provide a coupling between the impurity 
released spins. As we said, for larger doping, the background will be much more influenced by the 
impurities, which, besides creating states inside the gap, would also reduce the main gap, making 
a true antiferromagnetic long range order competitive with respect to the spin-Peierls phase, as 
discussed in the Introduction. 

Nonetheless, the neutron beam that probes AF ordering has a finite coherence length. Therefore 
the relevant question, even at low doping, is not whether a true long range order exists, but rather if 
clusters of size larger than that coherence length appear, and with which probability. For this reason, 
we want now to discuss AF ordering in a finite coherence volume, which, as we said, is set either by the 
size of the sample or by the finite coherence length of the experimental probe, e.g. the neutron beam. 
Therefore we vary the impurity concentration x at fixed size L. The Neel temperature is estimated 
by imposing a fixed value of the existence probability El(T) = 1/2. We are thus forced to analyze 
the low temperature behavior of the existence probability E^iT). This is shown on Fig. |7] for L = 40 
and impurity concentrations between 0.5% and 3%. With these parameters, the low temperature 
behavior of the existence probability is well fitted by the power- law 1 - E L {T) = T< x ). Fig. § shows 
the variations of the exponent a{x) versus x. In this figure, we have added the point {x = 0, a = 0), 
and, as it is visible, this extrapolation is compatible with the large x data, hence showing that zero 
temperature AF ordering occurs for an infinitesimal doping in a finite size. Imposing a fixed existence 
probability El(Tn) = 1/2, we obtain lnT/v oc —1/x, a form of the Neel temperature proposed from 



the low doping measurements in Ref. [14]. The calculated Neel temperature is shown on Fig. [|, 
together with the experimental points by Manabe et al. Compared to these experiments, we obtain 
a smaller value of the Neel temperature, still much above T*, in agreement with p4]j . 

In order to analyze the effect of increasing the coherence length, we now consider a system with 
a larger size L = 81. The low temperature behavior of the existence probability is shown on Fig. |9|, 
with doping concentrations from 0.25% to 1%, where the inclusion of a quadratic term was necessary 
in the fit. The Neel temperature obtained by imposing El(Tn) = 1/2 is shown on Fig. ||, together 
with the L = 40 result and the experimental data. As it is visible on this figure, in spite of the 
different behavior of the low temperature existence probability El(T) when the system size was 
increased from L = 40 to L = 81, the orders of magnitude of the Neel temperature T^{x) remain 
similar. 

The above analysis suggests that at low doping the model does not have a true long range order 
in the thermodynamic limit, but the probability to find an ordered cluster of any given size L is 
finite. This implies that the system will appear as if it were ordered to an external probe with a 
finite coherence length. This scenario gives a natural explanation of the fact that the measured Neel 
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Figure 7: Variations of 1 — E{T) versus the temperature T, E(T) being the existence probability for 
percolation of strongly coupled clusters in a system of size L = 40. The different curves correspond 
to: x = 0.03 and a = 1.25 (O); x = 0.025 and a = 0.98 (+); x = 0.02 and a = 0.72 (□); x = 0.015 
0.48 (x); x = 0.01 and a = 0.29 (A); x = 0.005 and a = 0.135 (*). The CuGe0 3 parameters 
44.7K,^ = 9and^ = 0.1^. 
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Figure 8: Variations of the exponent a involved in the fit 1 — E(T) = T a shown on Fig. |7[ The 
extrapolation to x = and a = has been added and is consistent with the data at higher concen- 
trations. 
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Figure 9: Variations of 1 — E(T) versus the temperature T, E(T) being the existence probability 
for percolation of strongly coupled clusters in a system of size L = 81. The fits are of the form 
log(l-£(T)) = /(log(T)), with x = 0.01: f{x) = -O.llx 2 + 0.18a; (x); x = 0.0075: f{x) = 
-0.092x 2 + 0.042x (□); x = 0.005: f(x) = -0.05x 2 (+); x = 0.0025: f{x) = -0.018x 2 (O). The 
CuGeOs parameters are A = 44. 7K, £ x = 9 and £ y = 0.1^ x . 



temperatures are always much larger than the typical energy scale T* of the doped CuGeC>3 , and of 
the evidence that the intensity of the AF peaks observed in neutron scattering is not saturated even 
at 1.4 K U for a Zn doped sample with x = 0.034. 

To obtain the dependence of the susceptibility upon temperature, we would need to improve our 
decimation scheme by renormalizing the exchanges, similarly to Ref. p^ ], and inserting a weak spin 
anisotropy, to break the spin symmetry. 

Let us summarize our decimation results by comparing them to the experiments by Manabe et 

al. 



14 , and to Tl 



Fit by Manabe et al: lnT N ~ 5.7 x 10 3 /x 
Decimation (L = 40,81): lnT^ ~ 10" 3 /x 
Behavior of T*: InT* ~ 10 _1 /x, 



(15) 
(16) 
(17) 



which illustrates both the fact that we can reproduce sucessfully the absence of critical concentration 
in our model, and the crucial role of disorder in enhancing the Neel temperature. 
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5 Conclusion and discussion 



Let us summarize the proposal we have presented for the onset of the Neel long range order in 
doped CuGeC>3. We first want to stress that, in our opinion, the real difficulty in the subject is 
not the coexistence of dimerization and antiferromagnetism, since the two order parameters are not 
incompatible. Instead, what is puzzling is the appearance of antiferromagnetism at such low impurity 
concentration as 1 x 10~ 3 , with a reasonably high transition temperature. We have shown that this 
behavior might depend on the peculiar properties of the disorder in this system, as well as in other 
spin gaped systems. 

We believe that the most important effect of disorder, at low doping, is to induce randomly 
distributed localized free moments. These moments are weakly coupled by an exchange constant 
exponentially decaying upon the distance, see Eq.(||), caused by the virtual polarization of the singlet 
background. We have shown that even a single chain does develop quasi long range Neel order, i.e. a 
power law decaying staggered correlation function, which, we believe, does imply a low temperature 
transition to an AF phase as soon as the interchain exchange is taken into account. Moreover, we 
argued that the exponential dependence on the distance of the exchange between those randomly 
distributed moments implies that rare events in which the impurity released spins are closer than 
the average distance (and therefore get strongly coupled), dominate the physics at low doping, in 
a Griffith-like scenario. This picture was fully confirmed by an investigation of the 2D disordered 
problem using a decimation method. The formation of these Griffith-like clusters, together with the 
finite staggered susceptibility of the singlet background, may explain why ordered domains larger 
that the neutron coherence length do appear at temperatures much larger than the temperature T* 
(the exchange between two moments at a distance equal to the average one). 

Some consequences can be drawn from this scenario: 

• A Neel transition is expected for any arbitrary small impurity concentration. 

• At relatively large doping concentrations, the mean field picture is expected to hold since the 
inhomogeneities of the AF order are weak. In this concentration range, the Neel temperature 
behaves like T/v — Jb^spx, where Jb is the exchange constant in the 6-direction, and £sp the 
correlation length in the pure CuGe03. 

• For low impurity concentrations, the Neel state becomes highly inhomogeneous and the dec- 
imation method is more suited than the mean field theory. Instead of a single transition 
temperature we find a quite broad distribution of T/v's. We have indications that a true long 
range order is never established, so that the low temperature phase might be rather interpreted 
as a Griffith- like antiferromagnet. We have shown that, for small coherence volumes, the Neel 
temperature behaves like In T/v ~ —1/x. By increasing the coherence volume, this behavior 
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does not rigorously hold. However, as it is visible on Fig. f|, the behavior of the Neel temper- 
ature for L = 81 can hardly be discriminated from the L = 40 calculation, and a fit of our 
results to a lnTjy ~ — 1/x behavior would be equally satisfactory. 

• A very long correlation length in the c-axis direction, extending well above ^sp, may be ob- 
served in a wide temperature range. This can be measured in NMR experiments at very low 
doping. For instance, the NMR linewidth in the ladder compound Sr(Cui_ x Zn x .)203 at a Zn 
concentration x = 0.25%, at which no Neel transition has been observed down to the lowest 
accessible temperature, is two orders of magnitude larger than the correlation length of the 
undoped compound p9| , which indicates that a similar behavior might well occur in the doped 
CuGe0 3 . 

We now compare our results with other theoretical proposals. We first compare our work to 
the one by Mostovoy et al. [27], and next to the work by Fukuyama et al. |35|], and Yoshioka and 



Suzumura |3q| . The technicalities involved in the discussion of these last two works are left in the 
Appendices [A| and [B[ and we focus here in the concluding section on more qualitative aspects. 

Mostovoy et al. |^7| have recently investigated the effects of interchain coupling at a mean field 
level in an antiferromagnetic Heisenberg chain with random dimerization, modeled by a (^-correlated 
in space gaussian-distribution peaked around a finite average value. In order to make contact with 
their work, we make the following comments: 

(i) The microscopic modeling of doping in a spin-Peierls system is quite different in our model 



and in Ref. [27|. As we showed in Ref. [|16|], our model is more related to a dimerization 
with a random telegraph noise of zero average distribution, with the low temperature physics 
dominated by the localized spins around the domain walls. 

(ii) The physics of the ID limit of both models is different. This can be seen by considering the 



RG flow [|n], 2C]. The ID limit of the model by Mostovoy et al. has been solved by Hyman et 
al. [ 28 1 , and flows to a fixed point with power-law Griffith singularities. By contrast, the ID 
limit of our model flows to the random singlet fixed point [^] with logarithmic singularities 



(iii) More importantly, the random dimerized Heisenberg chain, which is the starting point of the 



description in Ref. |27J, does not develop quasi long range AF fluctuations, as it is the case 
in our proposal. This allows to understand why a finite critical concentration of impurities is 
required in Ref. [£7|] to generate an AF-dimerized phase. As we show in the present article, 
this is not the case in our model, and it seems that this is not the case in experiments either 

0- 

A proposal similar to the one by Mostovoy et al. was previously put forward by Fukuyama, Tani- 
moto and Saito (FTS) |35|], who also worked out a model showing a coexistence of antiferromagnetism 
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and dimerization, as well as the onset of antiferromagnetism for arbitrary small doping. As we dis- 
cussed, we have started our analysis by assuming that the main effect of the impurities is the release of 
free spins out of the singlet background, which is in agreement with the magnetic susceptibility data 
on the doped compounds, showing an almost Curie like component at low temperature which seems 
to scale linearly with doping (see for instance Ref. B and [pf]). Instead, FTS assume, as Mostovoy 
et al., that the main effect of disorder is the local reduction of the dimerization 5i m p < 5 close to 
each impurity, hence of the spin gap. For instance, they can fit the experimental reduction of the 
scattering intensity from the dimerized lattice in the 0.7% Si doped compound, by taking <5j m „ ~ 0.25. 
This leads to a rough estimate of the local reduction of the spin gap A imp ~ (0.2) 2 / 3 A ~ 0.9A. Ob- 
viously, for smaller impurity concentrations, one should expect an even smaller reduction of the spin 
gap. However, to explain within this model the almost Curie like behavior of the low temperature 
susceptibility at low doping, one should assume a much stronger suppression of the spin gap down 
to ~ lOmK close to each impurity, not just a small reduction as the fit to the neutron scattering 
data would imply. This is one reason why we believe that our approach is more relevant. The other 
is that we do not believe that a local reduction of the dimerization might explain the appearance of 
antiferromagnetism at such low concentration as 0.1% of Zn doping. 

In fact, a more realistic modeling of disorder in CuGe03 should include both disorder effects, 
the local reduction of dimerization as well as the appearance of free spins. Essentially, the main 
difference between our approach and the FTS one is in the different emphasis placed on the two 
effects. We believe that the main role at low doping in establishing the Neel phase is played by the 
free moments induced by disorder, while FTS seem to favor the other mechanism of a local reduction 
of the dimerization. 



Apart from this comment on the grounds of the model proposed in Ref. [35], we also have some 
more technical criticisms. To analyze their model, FTS use the so-called Self Consistent Harmonic 
Approximation (SCHA) applied to the bosonized version of a dimerized Heisenberg chain. This 
approximation is a mean field theory, or saddle point approximation, plus random phase fluctuations. 
By this technique, they show the existence of a mean field solution having both a finite dimerization 
parameter and a finite Neel order parameter. Essentially, close to the impurity, where the dimerization 
parameter is assumed to be reduced, a staggered magnetic moment in the z-direction arises. We 
believe that this approach is not fully consistent, as we explain in detail in Appendix A. 



An alternative variational approach has recently been developed by Yoshioka and Suzumura [36] 
for the Zn doping case. They have rigorously shown that in bosonization the spin-lattice coupling 
does change sign crossing a non magnetic impurity, as it was proposed in Ref. [p^| . With this 
disorder modeling they have proved the existence of a variational action which does again describe 
a state where the dimerization gets reduced around an impurity and at the same time a static 
staggered moment appears in the z-direction. We believe that this variational approach gives a 
reasonable description of the Neel ordered phase, very close to what we propose, but does not allow 
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a description of how the AF ordering occurs, which we did in the present work. We postpone to 
Appendix |b] a detailed discussion of this point, and we end-up with two simple comments on these 
semiclassical analyzes. 

The dimerized phase breaks translational invariance like a Neel phase does, which, in addition, 
also breaks spin symmetry. If one analyzes the role of defects in a dimerized state within a scheme 
which breaks SU(2) symmetry (what is implicitly done in those semiclassical approaches, as shown 
in Appendix [B]), it is not a surprise that a finite staggered magnetization comes out of the calculation 
(without spin anisotropy or an applied staggered magnetic field), since the dimerized state already 
breaks translational symmetry. 

Finally, one can work out a model in the spirit of the FTS scenario. This can be done by 
considering a dimerization distributed as a telegraph noise, and taking a value 5 in the dimerized 
background (large segments), and < 5i mp < 5 close to the impurities (short segments). The 
resulting disordered model can be exactly solved in the XX limit along the lines in Ref. |57|. This 
would result in a model with power-law Griffith singularities, equivalent to the model by Mostovoy 
et al. and inequivalent to our model where the Dyson singularities are power-law. This model 
would not flow to the random singlet fixed point [20 J, would have much weaker AF fluctuations than 
our model, and, as in the case of the model by Mostovoy et al, a critical concentration would be 
expected for the appearance of AF order. 



6 Acknowledgments 

We acknowledge useful discussions with J.C. Angles d'Auriac, C. Berthier, D. Fisher, T. Giamarchi, 
B. Grenier, M. Horvatic, J. C. Lasjaunias, J. Lorenzo, P. Monceau, A. Nersesyan, Y. Okabe, C. 
Paulsen, J. P. Pouget, G. Remenyi and M. Saint-Paul at various stages of this project. The authors 
also thank K. Uchinokura for sending them their susceptibility data. This work has beeen partly 
supported by INFM, Research Project HTSC. The decimation calculations have been performed on 
the CRAY T3E supercomputer of the Centre de Calcul Vectoriel Grenoblois of the Commissariat a 
PEnergie Atomique. 



A Self Consistent Harmonic Approximation 



In this Appendix, we discuss the Self Consistent Harmonic Approximation (SCHA) approach used 
by Fukuyama, Tanimoto and Saito (FTS) in Ref. [pq |. We do not want to enter in the details of 
the transformation which allows to map the spin-Peierls problem onto a bosonic one, since they can 



be found in Ref. 3S], and references therein. We assume that this transformation is valid, and, 
following |j3q1, that the resulting bosonic Hamiltonian is 



H 



dx 



U(x) 2 + (d(j)(x)) 2 - 2gu(x) sin (V^ttK^x)) + 



(18) 
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where 4>(x) and n(x) are bosonic conjugate fields. The classical variable u(x) represents the dimer- 
ization lattice distortion, and the sin(. . .) is the spin dimerization order parameter. Notice that the 
cos(. . .) would represent the z-component of the staggered magnetization. We then write 4>(x) = 
<t>c(%) + fifa), where 4> c (%) is a classical variable, and 4>(x) is the quantum fluctuation part, which is 
assumed to have a zero average value. We want to find the minimum of the total energy assuming a 
bosonic state |0) which is the ground state of the following Hamiltonian 

1 



dx 



U(x) 2 + (d(P(x)y + m(x)(p (x) 



(19) 



with m(x) > 0. Therefore we have to find the appropriate values of <f> c (x), u(x) and m(x) which min- 
imize the average value of H over |0). This is the so-called Self Consistent Harmonic Approximation 
(SCHA), which is a well defined variational technique. Since ( |l9| ) is quadratic, then (0|^(x)|0) = 0, 
and 

(0| sin (yjAirK (<f> c (x) + </>(x))) |0) = D(x) sin (V^ttK^x) 



where D(x) = e 2llK (4>( x ) 2 ) _ jf we define T{x) = (II 2 (x) + (d(j)(x)) 2 ), the self-consistent equations are 



-d 2 <ft c (x) — gV AttKu(x)D(x) cos (^/AnK(p c (x] 



KU(X 



gD{x) sin (V 4nK(j) c (x) 



5m(x) 



(20) 
(21) 

(22) 



where we used 



We start by solving 



6D(y) 



-2nKD(y) 



5{4>{yf 



5m(x) yt " 5m{x) 

). We write the bosonic fields as 

t 



where if- a (x) = ip* a (x), e a = e_ a , 

Y,vl( x )<Pa(y) =S(x-y), 



(23) 
(24) 



ti 



and [b a , bp 



d a} 0, all other commutators being zero. We want the Hamiltonian (19) to be diagonal 
in terms of these creation and annihilation operators: 



A straightforward calculation gives 



d 2 ip a (x) + m(x)(p a (x) = e 2 a ip a (x) 



(25) 
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which is the eigenvalue equation for (p a . On the bosonic vacuum of the operators b a , we find that 



e« V* a {y)v a {y) + —d<p* a (y)d<p a (y) 



J 'dyT(y) = i]T J r dy 

which, through Eq,(p5|), becomes 

J dyT(y) = J dy (2e 2 a - m{y)} ip* a {y)ip a {y) 



E 



E €a ~ / d y m {y) l \4> 2 {y)), 



where the last identity follows from (f23|). Now, suppose that m(x) — > m(x) + <5m(x). It is easy to 
get the first order correction to the eigenvalues through (p5|): 

1 



5>* = E 



rfy <p* a (y) $ m (y) faiy) 



2e,__ 

dySm(y)((j) 2 (y)) 



{(j) 2 {x)) + J dym{y) 



Through the above equations we find that 

J dm{x) 
which, inserted in Eq.(|2"2|) gives 

dy —m(y) + AirKgu(y) sin (V £irK4> c (y)j D(y) 
This equation is solved by imposing 



H<P 2 (y)) 

5m(x) 



H4 2 (y)) 

5m(x) 



0. 



m{x) = 4irKgu(x) sin I v 47rK(j) c (x) J e 



-2irK{cf>(x) 2 ) 



(26) 



which gives the self-consistent equation for the mass. 

As a simple application, let us consider the homogeneous case m(x) = m, u(x) = u, and 4> c (x) 
(p c for any x. From Eq. ((20) we get 



y/4xK<j> B (x) = 7r(n+ -), 



(27) 



which choice is in fact the only compatible with SU(2) symmetry. We choose n = 1. On the other 
hand Eq.(|25|) is solved by plane waves 



ip q (x) 



1 



7F 



-iqx 
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with energy e q = \/ q 2 + m. Therefore 



IX! \2\ 1 1 1 i ( m 



m 



where mo is a high-energy cut-off, and ( |26[) becomes 

/ m \ K/2 

m = AirKgu ( I 

\m J 

The solution is m ~ (gu) 2 ^ 2 ~ K \ Since the approach is valid only if, for a small coupling g, a small 
mass m is obtained, the above equation implies that a mass is generated if K < 2. This is the famous 
result for the sine-Gordon models. If we consider this model as describing a spin-Peierls transition, 
we know that K = 1 corresponds to the XX-limit and K = 1/2 to the SU(2) point. In the latter 

case, 

/ m \ 1/4 
m = 2irgu — 
\m J 

Eq.(pl|) is solved by 

g f m\ l / A 
K \m J 

finally leading to 



m= 4 4 l - (28) 



Now let us consider the inhomogeneous case. For instance, following FTS [35], we assume that 
■u(O) = u(l) = Ui. Through (p6|) the self consistent equations can be written as 



-d 2 <p c {x) - gV4nKu(x)D(x) cos (y^rK A c (xfj = 0; (29) 
ku(x) = gD(x) sin (V4irK(j) c (x)) ; (30) 



m(x) = 4irKgu(x) sin yj AitK<t> c {x))D{x). (31) 

Solving for u(x), we get 

and the coupled equations for m(x) and (j) c (x) read 

^^D 2 {x) sin 2 (V4^0 c (x)) , (32) 



m{x) 

K 



(33) 



- d 2 M*) ~ g2Vi ? Z)2(x) sin (2v^K0 c (x)) = 0. 

We do not know whether these equations have an unique solution. FTS implicitly assume that D{x) 
as well as m(x) are independent of x. Then Eq.(|30|) at x = 0, 1 fixes the value of C (O) and (j) c (l), 
which are used as boundary conditions to solve Eq. (|33|) . The solution is a Jacobi elliptic function. 
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Between x = and x = I, the solution approaches the homogeneous one (|27|), but around these 
points it is different. As a result, sin 2 (2\J Ait Kcj) c (x)j differs from one close to x = and x = I, 
leading to a finite value of cos 2 \ 2y/ 4irK(f) c (x)j , hence to a finite staggered magnetization. This 
result would imply that disorder in the lattice distortion leads automatically to a long range Neel 
phase, but necessarily with a staggered magnetization along z. In fact, the staggered magnetization 
perpendicular to z is expressed in terms of the conjugate momentum H(x), which is uncertain if 
<f>c(x) gets an average value. This result is puzzling, since the starting Hamiltonian is spin isotropic. 
Moreover, the assumptions D{x) and m(x) ^-independent have to be aposteriori verified. FTS do 



not make this check, but clearly, by inspection of Eq.(32), one realizes that the assumption is not 
correct. Therefore the calculations of Ref. are not self-consistent. 

Alternatively, one can look for a solution which does not break SU(2) symmetry. In this case we 
have to assume a constant (fr c (x) = 7r /2, compatible with (E7l). Therefore one needs to solve 



ku(x) = gD{x)\ (34) 
m(x) = 2irgu(x)D(x), (35) 

or, equivalently, 

2ir 

m{x) = — g 2 D 2 (x) = 2ttku 2 (x), 

with the boundary conditions m(0) = m{l) = 2ttku 2 . If this equation were solvable, we would get 
an alternative solution which keeps the SU(2) symmetry. Further work to establish which is the 
best self-consistent solution is needed, but we can safely state that the assumption of homogeneous 



quantum fluctuations made in Ref. [35| is not self-consistent 



B Some rigorous results in the presence of solitons 

In order to better clarify the limits of the SCHA, we consider a simple case which is related to the 
work by Yoshioka and Suzumura [36]. These authors have proved by bosonization that the coupling 
to the lattice changes sign crossing a non magnetic ions. A similar conclusion was reached in Ref. 
[ |16| , by means of a simple mapping to a squeezed chain where the non magnetic sites are eliminated. 
This mapping is explained in Figjl). The squeezed chain is obtained by eliminating the non magnetic 
site similarly to what is done in the U —* oo Hubbard model. The difference here is that the hole 
does not move. Next, the weak link across the impurity is approximated as the weak bond of the 
dimerized pattern, i.e. J (1 — 5), which is supposed not to change qualitatively the low energy physics. 
Hence, the effective model describing the squeezed chain consists of a dimerized Heisenberg chain 
with a domain wall. The model in the presence of a single domain wall can also be analyzed along 
the same lines outlined in Ref. |36|. The result would be again that a variational solution can be 
obtained with a space dependent classical phase 4> c (x), but a space independent mass for the bosons 
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describing the quantum fluctuations. The same would hold also in the XX-limit, which can be also 
easily solved by diagonalizing the spinless fermion Hamiltonian which is obtained through a Jordan- 
Wigner transformation. In the continuum limit and after linearizing the band around the Fermi 
momenta kp = rbrr/2, one has to solve the following equations 

d 

- i-r^XeR(x) ~ iA(x)XeL(x) = €XeR.{x) , (36) 

d 

i-Q^XeL{x) +iA(x)xelt(x) = £XeL(x), (37) 
where the Fermi field is expressed in terms of the x functions through 
*(x) = [e ikFX Xen(x) + e- lkFX X eL(x)] c e ^Y,< 



[X )C f 



c e being the annihilation operator of an eigenstate of energy e. The dimerization parameter is equal 
to A for x < 0, and to —A at x > 0, thus having a domain wall at x = 0. Eqs.(|36"|) and ( |3~7| ) can 



be solved and one finds, in addition to scattering solutions with energy e = ±V& 2 + A 2 , k being the 
wavevector, a localized zero energy solution with wavefunction 

Mx) = V2Ae- lxlA sm(k F x). 

This is a soliton solution for a single domain wall. 

The classical phase which appears in the bosonized version of the model can be related to the 
phase of the 2kp oscillating part of the spinless fermion density. Namely, we expect that 

P2k F (x) = {^ j (x)^(x)) 2 K F ~ cos (2k F x + (j) c (x)) . (38) 

The dimer order corresponds to a bond ordered Charge Density Wave (CDW), i.e. <f> c (x) = (2n + 
l)7r/2, while the Neel phase to a site ordered CDW, i.e. <fi c {x) = nir. In the absence of the domain 
wall, the phase is <j) c {x) = (2n + l)vr/2 at any x. In addition, in the presence of the domain wall, two 
other terms arise. One is 

pxix) = (n - Ae- 2 ^ A (1 - cos(2k F x)) , (39) 

where no is the occupation probability of the soliton state. In fact, the term proportional to no is 
the soliton wave function density probability, while the other one comes from the negative energy 
scattering solutions. We notice that p\ oscillates with a phase which has two possible values cp c = 0, it 
depending on no, which are compatible with the two possible Neel ordered phases. The other term 
which is generated by the domain wall decays like 

P2 (x) ~ ^ 6-^, (40) 
v Ax 
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and generates the polarization of the background responsible for the coupling among consecutive 
solitons, which we have considered throughout this work. This term oscillates with a phase compatible 
with a dimer order. Were pi(x) be present, we would indeed find a coexistence of dimer order and 
Neel one, the latter mainly localized close to the domain wall. 

We have now all the elements to describe what happens in this solvable model for the classical 



phase, and compare with the variational approach of Ref. 36]. We start reminding that after 
the Jordan- Wigner transformation, a site occupied by a spinless fermion corresponds to a the spin 
having a z-component S z = 1/2, while an empty site to S z = —1/2. Moreover, a chain with periodic 
boundary conditions is compatible with the presence of a single domain wall only if it has an odd 
number of sites. Therefore the ground state has a total S = 1/2. In our XX-limit we can not 
rigorously discuss the SU(2) symmetry, but at least the following discussion gives some hints how the 
implementation of such symmetry gives a classical phase <p c (x) = vr/2 independent of the position. 
In fact, if we take a ground state which does not break the symmetry, we have to assume that the 
chemical potential crosses exactly the soliton energy, so that there is an equal probability to have 
S z = ±1/2. In other words this implies that no = 1/2, leading to a vanishing pi(x), and leaving only 
the P2(x) term which has indeed a phase fixed at 4> c {x) = ir/2. Let us suppose a domain wall and an 
anti-domain wall at a distance much larger that 1/A. We can now consider a periodic chain with an 
even number of sites. The ground state will be a singlet. Around each domain wall there will be a 
localized soliton state at zero energy. The two S = 1/2 solitons will be coupled by the polarization 
of the background, generating a singlet S = as well as a triplet 5 = 1 state. The singlet state 
has energy —t (0 < t <C A), the S z = component of the triplet +t, and the S z = ±1 have energy 
zero, where t decays exponentially with the distance among the two domain walls. The fact that 
the triplet is not degenerate is a consequence of the XX-limit. In this case, however, we are lucky 
since the lowest energy state, being a singlet, is SU(2) invariant. This is in essence why the spin 
anisotropy with a random distribution of domain walls does not play a very relevant role for what it 
regards the ground state properties. Within the lowest energy singlet state, each soliton state has an 
equal probability of being empty or occupied, leading again to a vanishing pi(x) Eq.(p9|). However, 
one can view the singlet state as oscillating between two configurations each having p\{x) ^ 0, hence 
a classical phase <j) c = 0, it compatible with a Neel order. In a rough picture, the system is locally, 
in time, in one of the two possible Neel ordered states, and oscillates between them over a very 
long time scale ~ 1/t. Moreover, in spite of the exponential decay of the polarization over a length 
S.SP ~ V(2A), which is the correlation length of the dimerized state, the two solitons, which are at a 
distance much larger that £sp, oscillate coherently. This is on a very small scale our picture of what 
happens in the doped CuGeC>3. 



We think that the variational approach of Ref. [36] is able to describe only the states where 
one soliton is occupied and the other is empty, or viceversa, and not the superposition of the two, 
which is beyond a mean field like approach. Therefore the method implicitly breaks SU(2) symmetry, 
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which obviously leads to a finite staggered magnetization pi(x), since already the dimerized state 
breaks the translational invariance. For this reason, we believe that this approach should give a 
reasonable description of the Neel ordered phase, but is unable to describe the way in which the 
magnetic ordering is established, which is what we have done in this work. 
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